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ABSTRACT 


The structure of reel numbers is too rich to admit 
their proper representation even on idealized, computer such 
as an Unlimited Register Machine, It follows that if a 
physical problem is formulated in terms of real analysis, 
then there would be on unbridgable gap between what 
analysis would show about the nature of the physical problem 
and what a computer simulation T would yield about it. An 
attractive alternative is to switch over from real numbers 
to what are called computable numbers. The analysis based 
on computable numbers centers .around the fact that a number 
exists if and only if it can be computed. 

In this thesis an introductory account of computable 
numbers and computable analysis is presented for the readers 
concerned with digital signal processing. This is done with 
a view to see how far it would be feasible and worthwhile to 
reformulate concepts and algorithms of digital signal process 
ing in terms of computable analysis. 
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Section 1 


Introduction^ 


1.1 Mot ive bion. - Real analysis v/s computable analysis : 

In trying to introduce the notion «f compu table 
numbers and explain what purpose is served by tnem, it is 
perhaps best to start with the example given by Aberth [ l] . 


Let a function be defined as follows s 


f.(x) =0 if x is a rational nurciber 
f(x) =1 if x is an irrational number 

Now let us see, what the value f(x) assumes, when x equals 
the constont 7, where 7 is the Eulerfe constant and defined as 


7 
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At the present time it is not hnorm whether V rational or 
irrational, so we cannot assign a value to £ (7) . 


In real analysis this difficulty is commonly exp- 
lained away by saying that while the function does exist 
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Our not being able to compute its value simply reflects the 
imperfect state of our knowledge. 


But then what does it really mean to say that some- 
thing does exist but we cannot exactly ascertain what it is. 
One could take the view that a number exists only to the extent 
that we can compute it. It is this view which has led to the 
notion of computable numbers/ functions or sequence. Detailed 
description about computable analysis and computable number 
will be given m section 4 f but here ua shall see what we mean 
by a compu table number m very broad terms. 

Let us begin by considering a simple example illus- 
trating the fact that if we were to consider computation 
in traditional manner we would perhaps look at series expansion 


x 

e 


' 2 3 

i x , x 

1 + * + H + "3 


* • *- 


and say that determining its value is only a partial sum of 
the series. But then this heaves en entirely open question 
about errors. If we wish to incorporate here the fact that 
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its computation be possible using a finite programme in a 
finite number of steps, it is not enough to just give the 
series expansion. What we need in addition is to specify 
on error bound such that it can be ascertained within this 
error bound by an algorithm. We will use X(s) to denote the 
computable numbers approximation of X which has an error of 
t As we shall see, is a number that can be determined 

<r«. <3 

by computation On an idealized computer and the set of values 
X(s) for different values of ‘s’ as the outputs of the same 
computer program producing these values for different values 
of the parameter s. 

As we will see later also in detail a computable 
number can be taken in general to be a programme or algorithm 
of the kind that can compute a function within the required 
error bound and terminate execution after finite computations » 

Now to illustrate the difference between computable 
analysis and real analysis, we will consider two elementary 
results of analysis ’• 

1. The first one is the Rolls’ s theorem. It states that if 
f(x) is a continuous function on a closed interval a,b 
with f (a) and f (b) of opposite signs, there is at least 
one real number C in a,b such that f(c) = 0. Now if we 
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switch over to computable numoers then it can be shown that 
this theorem docs not hold. More precisely it can be shown 
that if f(x) is a "computable function on the closed interval 
a/b with £(a) and f(b) of opposite signs, then there may not 
be a computable number C in a,b such that £(c) =0 [8] . 

2. The second result is that function f(x) continuous on a 
closed interval a,b assumes its supremum at a point in 

4 

this interval. The corresponding computable analysis 
result is that a corrputable function f(x) continuous on a 
closed interval a,b may not assume its supremum at a 
point in this interval. 

Thus there are areas of differences between compu- 
table and real analysis. 

1 . 2 Dig_ital_ Sig nal Proce ssing 

Signal processing is a field in which we use programs 
to process the raw data or for designing digital signal pro- 
cessing hardware systems. Theoretical results in digital 
signal processing are primarily based on real analysis. But 
when models based on real analysis have to be implemented on 
computing machines, there are serious problems owing to limi- 
tations of machines in representing real numbers. Thus there 
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a ground to think that if one uses computable analysis/ it can 
help evolve a more realistic theoretical framework for digical 
signal processing because computable numbers admit exact repre- 
sentation on idealized computer. 

The aim of tins thesis is to bring together tnose 
basic concepts and results of computer analysis m an intro- 
ductory fashion that are of potential use in digital signal 
processing. The original plan wes to carry out reformulation 
of the theory of digital filters in terms of comp u table numbers/ 
however/ as work progressed, it was felt that the task was 
not as straightforward. There were independent areas of work 
connected with computable analysis that needed to be studied 
and sorted out. The iwork reported here is the result of such 
a study. 

3.3 Ove rvi ew o f S ections 

As the aims of thesis are to represent computable 
number concept and connect Digital Signal Processing with it, 
sections 2 to 4 have discussions centred around computable 
number/ computable analysis and computable arithmetic. Section 
2 e:xplains the need of computable numbers in place of real 
numbers for computational work. 
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Section 3 is concerned with general discussion on 
computable analysis. Here the so called most accurate possible 
(MAP) real analysis is compared wich computable analysis. 

A orief historical account is also given. 

Section 4 is about computable number arithematic 
operations and relations with computable numbers are explained 
and the fact that computable numbers form a field is discussed 
in detail. 


Section 5 and 6 explain the use of computable numbers 
in Digital Signal Processing. Section 5 contains an account of 
the modification required to make computable analysis a practi- 
cal analysis and atterrpt is also made to show the possible 
algorithmic approaches to implement the required modifications. 

Section 6 deals with the digital signal processing. 

An effort has been made to reformulate the basic properties 
of a DSP system in terms Of computable numbers. The Dead band 
effect and problem of sensitivity are considered. 

Section 7 contains, in addition to concluding reamarks, 
certain suggestions for further work. 



Section 2 


Ap p r oqr ia-toncss of Computable Numbers 

2 .1 T he Ex cessive Richness of Real Numbe rs 

In the traditional hierarchy of numbers we start with 
natural numbers followed by integers, ratxonals and ultimately 
arrive at the real numbers. In the real number system we 
reach a state of prefection which cannot be excelled so to say. 
Among the various mathematical entities that we hnow so far, 
real numbers have the richest structure possible. Further, 
for purposes of modelling and analysing physical systems, the 
real number system furnishes a foundation which is well tested 
and of long standing. 

But this enormous richness of the structure of the 
real number system is in some ways a hindrance to a realistic 
modelling of physical processes and systems. In the early days 
of scientific development when the scale of physical problems 
was not so large and when analytical models and solutions were 
considered adequate, this fact did not draw much attention. 

In the modern setting, however, when we wish to 
model, represent and analyse very large systems purely by 
computational means, the use of the real number system and 
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mathematical analysis based on it cannot be accepted without 
reservations. There are two very basic issues involved here, 

Firstly the concept of real number as we conunonly 
understand it todays do^s not contain within it any reference 
to computation. We are concerned -with the fact that it exists/ 
but we are not concerned whether there is a particular algorithm 
or construction by which its value can be determined with a 
given accuracy and with a finite amount of computational effort. 
Yet when our objective is to simulate mathematical models of 
physical processes computationally/ then one cannot meaningfully 
conceive of a number except in terras of an algorithm that 
produces it. There is thus the basic limitation with real num- 
bers that there may be those that cannot be defined in terms of 
effective computational algorithms. In case it is SO/ then 
there would remain on ubridgable gap between the mathematical 
models and the computational simulations for tnem. 

Secondly/ large scale computations require the use of 
computers which when reduced to their essentials/ are finite 
state machines constrained by limitations of memory and speed 
in a manner that does not permit an exact representation of 
computations with reals. 

Traditionally/ these limitations are handled by using 
what is known as most accurate possible OXIAP) representation of 
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real number [ 7 ] . In this representation, a real numb or 

is represented as accurately as possible on a given computer 
and the computations are carried out maintaining the highest 
accuracy possible. It is hoped that end result would tnis way 
automatically have the necessary accuracies. A question that 
arises is the following s 

Instead of 'fitting' the real number system on to 
a computing machine can we not work backwards. That is, 
starting with the structure of commutations as they are carried 
out on a computer, can we choose a number system that is exactly 
implcmentable atleast on an 'Idealized' machine and is also 
close enough to the real number system. 

In regard to this question, we find from the literature 
[ 1 ] that the computable number system is one such choice. 

In this number system every reel number is approximated by 
computabl e values lying within the errorbound specified in 
advance. This system acknowledges the limitations of the com- 
puter in representation of real numbers and treats real numbers 
as knowable only to within certain finite errors by means of 
computational algorithms. 

In order to fully appreciate the significant of 
computable numbers, it is perhaps essential first to have an 
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idea of Idealized computing machine models and their limitations/ 
and to see how computable numbers senre as a matching device 
between computers and the real number system. 

2.2 Idea lized Computer Ilodels 

In the literature/ ideal computer models range from 
Turing machines to some higher level models as Unlimited Register 
Machine [2] . All these modcjls rely on the availability of 

infinite memory in one wav or the other. On the ground of the 
thesis [2] that all these de different models are computationally 
equivalent/ we propose to confine ourselves to the URM model. 


The URM model consists of an infinite number of regis- 
ters labelled Rl, R2, R3, ... etc #/ each of which contains a 
natural number. This loolcs as shown in fig. (2.1 ) 


Rl R2 R3 R4 R5 ... 



Fig. 2.1 Representation of a URM 

For computations on URM, one has to write a program 
containing instructions from a set of standard instructions 
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which a URM Can execute. The standard instructions are s 


1. Z ero Instruction s By executing chis instruction URM puts 
‘O' in the nth register. It is donated by Z(n) for n = 1,2,3/... 


etc. 


Example Z (4) 






2. Success or . In^-s trucj: l p n : This instruction increases the 

number contained in by 1 and is represented as S(n) for 
n ~~ I, 2, 3 ... etc* 


Example t S(3) 



3. Transfer In s jtruction : This instruction, represented as 

T(m,n), replaces the contents of R n by that of R m , for m, n = 

1,2,3, ... etc. 

Example : T(3,4) j~"2 ' 5 ~ L ^ I 4 [ 3-L 

4. Jump Inst ructi on s This instruction branches off execution to 
the p th instruction if contents of R m equal the contents of & n 
and is represented as J(rn^n,p) where m,n,p = 1/2,3, ... etc. 

Example J(4,6,20) - if contents of R 4 = contents of R g 
then go to instruction no. 20, 
otherwise go to next instruction. 
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An example of a URM machine programme for computing (2) is 
the following i 


1 Z (4) 

2 Z(5) 

3 Z ( 3) 

4 S ( 3) 

5 S (4) 

6 j( 3/ 2,8) 

7 j(l,l,4) 

8 S (5) 

9 17(5,1,11) 
10 J (1,1, 3) 


2 

nn 
u n _ 

3 

LlJ 

0 

Rl 

R2 

R3 

R4 

R5 


This example shows how a function can be represented 
by a program of finite instructions. 

At this point a word about the provision of unlimited I 
memory in the URM model would be appropriate, if we assume that 
one instruction can be put in one memory register then a program 
in general can be written in a finiue number of memory locations. ! 
But there is no upper limit to the number of locations that 
may be utilized in a particular program. The unlimited register 
model is assumed to ensure that any finite step program, no matter; 
how 1 arge, can be accomodated in the memory . 
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Another assumption made for the URM model is that the 
word length of the computer be unlimited* This again is to 
ensure that integers, no matter how large, can be -stored. 

After examining this idealized computer model, we now turn 
to computable number and real number representations and 
their suitability. 

2.3 Representing Numb ers on C om puters 

About the idealized computer model, a natural ques- 
tion to ask is whether an exact representation of numbers 
system is possible on it. If it so turns out that a given ^ 
number system cannot be exactly represented even on an idea- 
lized model then the question of being able to represent it 
on a practical machine does not arise at all. The point that 
goes strongly m favour of computable numbers is that while 
not all real numbers have an exact representation on a URM or 
any equivalent idealized machine, all computable number do 
have such representations by definition. Computable numbers, 
as we will see in Section 4, by their very definition consist 
of programmes on a URM. Two specific properties of computable 
numbers that are of fundamental significance in this context 
are the following* 

1* A conputabl e number consists of an effective algorithm 
that can be realized as a programme on a URM. 
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2. Computable numbers like the real numbers have the structure 
of a field/ a fact that is known as Rice theorem. 

In the light of these two facts it is reasonable 
to say that computable numbers provide theor stical basis for 
computer based system studies that is more sound than the 

traditional one based on real numbers. 

2 .4 C omput able H umber <- R ecursive or Co nstructive : 

Before we close this general discussion about the 
appropriateness of computable numbers/ it is wortn while examin 
ing their nature in terms of the notions of computability and 
constructibility. This discussion is based on the idea that 
all recursive numbers are not constructive [3]- 

Here ‘ recursive* means while constructing the number/ 
a computer computes its first decimal digit then the next and 
so on. The programme also has some condition to terminate the 
execution . But it is not certain that the condition will ever 
be met. The term 1 Constructive 1 means that the condition to 
terminate the execution of programme will be met after a finite 
number of computations. Myhill has given examples of two recur 
sive number sets that are not constructive and one that is also 


cons true t ive 
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The first one is 'R ^ 1 t the set of decimal expansion 
of real numbers. Such the decimal expanded form terminates at 
the n t * 1 decimal where ‘n‘ is the ‘critical number of tt*. Here 
the critical number of n nreans the number places after which 
the string 7777 appears in the decimal expansion of it. So far 
it is not known what this number is. 

To give a specific example/ let us construct a number/ 

X/ from the set R^/ as follows/ we put x = 0.3333 ... t and 
continue repeating 3 till we reach an odd place 2n-.-l such that 
2n-KL is the* critical number of it 1 , then put 4 at 2n4-l° place 
and repeat the process. 

The second one is R-]_/ "t-h® s®" 1 - a H * locatable 1 real 
numbers. A number of this set can be constructed as follows. 

Consider a real number x and an integer B serving as 
the upperbound limit. Then compare x with -B, -3+1 , ... 0/ ... 
B-l/ B. Suppose we find by the comparison a number k such that 
B+k is less than x and B+k+1 is greater than x. Let x be a number 
whose integer is B+k. For the decimal part, coitpare x with 

2^ for n ~ 1/2/3 ... and find s such that B+k-h'^j~Q , x < 

B+k+ The 3®®!™® 1 pl® ce o£ x is occupied by s. Next 

we carry out similar comparison with B+k+ r § in place of 3+k. 
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This gives the second place of decimal. Repealing this defines 
the infinite decimal expansion of x. 

The third one is R^ the set of decimally approximate 
real numbers and whose elements are such that for any given 
rational e > 0 we can find a finite decimal expansion d such 
that x - d < e then d ^ R . 


Myhill has given theorems stating that and are 
not closed under addition and R^ forms a field. Thus decimally 
approximated real numbers are likely candidates for a construc- 
tive analysis but not or R d . The set R^, and R^ are all 
recursive, but R da is the only one that is constructive as well. 
Based on this observation of Kyhills, we can say that in working 

with the set R da we are dealing with computable processes, where 

1 

the error e of R. corresponds to the error - of a computable 

C13L S 

number X(s). On further comparison we find that the field R da 
is the same as the field of computable numbers. We thus conclude 
following Myhill that computable numbers are recursive as well, 
as constructive in the sense defined above. No w we can enume- 
rate a few examples of computable numbers which will strengthen 
the discussion. 
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2.5 A F ew 5 xarnp l.es of Computab le Number s 

To fix the idea discussed/ so far we list here 
-fewo examples of computable numbers. We will see from 
these examples that while the set of computable numbers is 
smaller than the set of reals/ the computable number set includes 
all numbers of practical interest. 

1, All rationals are computable numbers. 

2. Two computable numbers when added/ subtracted, multiplied 
or divided again yield a computable number. 



Section 3 


.About Congou table_ ^ analy sis 

3.1.1 Int roductions 

Having seen broadly what computable numbers are and 
what their place is in relation to the classical number systems* 
the next natural step is to examine what is meant by computable 
analysis . 


As real analysis m its classical form does not take 
into account the limitations m computations on a machine, we 
have two alternative courses as described m the previous 
section, which we can summarise as follows s 

1. Real analysis be used but real number represents- r , 
tion on machine m the MAP form on that machine. 

2. We resort to a representation of numbers which is 
matched to suit the limitations of the machine on 
one hand and has all the properties to admit a 
complete analysis on that number system analogous 
to real analysis. 

What we get persuing the first alternative is the so 
called Most Accurate Possible (MAP) real analysis, and the second 



alternative consists of computable analysis. Here vs will 
compare these two approaches and will give our arguments why 
computable analysis be favoured. 

We will also point out the interconnections of compu- 
table analysis with developments in the areas of con atun ict ive 
analysis and interval analysis. It is our thinking that studying 
these connections is important from the point of view of making 
computable analysis practicable. It will also help in parti- 
cular t in clarifying the fundamental fact underlying computable 
analysis that a number is a program. 

3.1.2 Dif f loul ti as with, MAP _real_ Ana lys is 

To be gin with, let us concentrate on what are the 
difficulties which make MAP real analysis inferior to computable 
analysis. MAP real analysis has the following four well known 
pitfall s. 

1. Real number is represented with some approximation 
because a machine cannot represent an infinite 
decimal real number value. So an error exists at 
each computation and this error part goes on buil- 
ding up as computations continue. 
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Let us say for instance we wish to compute 


,/ \ X , X 

(x) = x - + r-c 

1 ‘- 1 £ 


+ 


X 


2n-l 


]2n~l 


Let y be the lAP reprosen cation of x, such that 
y - e < x < y + s 

where e is the error bound, then in computation for 

F(x) in I LAP real analysis first term introduces the error + 

3 

Second term -pr has error bound 
ire 


\ 


x: 

U2_ U 


1 2 „ 2 3. 

t (3y s+ 3ye + e ) 


e sum(2) " £ ! + e 2 “ £ + I (3y2 £ + 3yt ^ + e3) 


This way we may compute the summation of n terms and. 
the corresponding error « sum ( n ) O' 0 ® 5 on increasing with n. 
Following the MAP approach there is no way for us to exert any 
meaningful control on this error. 

2 . We obtain two distinct results for the same function when 
it is computed on two different computing machines, e.g. , a 
hand claculator and a 20 decimal place precision computer. 
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3. Sxnce 11AP real analysis is not error conscious, certain 
results may be obtained 'which are totally unacceptable. 

1 


Let us compute 


5(1 -i- 10" 23 ) - 5 


using 


(a) ten decimal precision (b) 30 decimal precision 


Solving with ten decimal precision 


5(1 ; 10 - 5 


" 5 ( 0 . 999999 9 999 - 1)" 

1 _ 

5 x 10" 

Now calculation using thirty decimal precision 


1 1 .10 

r-Td - = - s x 10 


— (A) 


4. 


X 


5(1 + 10“ 23 ) - 5 


= + i x 10 23 


5 x 10' 


-2 3 


— (3) 


o 

We see that results (A) and (B) are too divergent. 

One refinement of MAP analysis is interval analysis, in 


which in addition to the MAP confutation an error analysis 
is simultaneously carried out. Even in this refined form, 
one has no way of controlling error to the required value. 
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All these four difficulties with the MAP real analysis, 
are a result of the fact that there are errors in corrou cation 
and nothing can be said in advance about their values. One has 
to just compute to the highest accuracy possible and hope for 
the best. If we wish to improve the situation then we have to 
look for some form of analysis which involves the concept of 
error bound in the number representation itself. Efforts in this 
direction leave given rise to what is known as cormutable analysis. 

3.1.3 Getting over jche J)i , f £ lcul _t oes 

As we shall see in section 4, a computable number is 
always thought of in terms of an error bound as a parameter. 

Thus we speak of a cornpuuable number X(s), which represents the 

. . 1 

real value X, where the value of X lies in the rarue X(s) - — 

s 

1 

to X(s) -l- r . If we think of numbers in this manner,, then 

« 

difficulties faced in MAP real analysis can be overcome one by 
one. 


Difficulty (1) with MAP real analysis is removed as 
the final error bound on the result is fixed (+ g- here) before 
we start execution of the program. So the error build up owing 
to large number of computations will not be more than the required, 
error bound. Thus if we wish to compute 

X(100) = Sin 100° 
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By the summation of series expansion of Sin 100°, 

is to foe continued upto the point that the truncation 

1 

error is less than or equal to Jqq* Computable number concept 
thus ensures a fixed error bound on the result. 

Difficulty (2) Is due to different precisions of 
different computing machines. But a computable number computed 
on two different machines will have the same value because both 
the computable numbers will have equal error approximation to 
a real number. 

y (100) = Sm 20° * log 5.31 

represented computable analysis version of the real number 
y = Sin 20° * log 5.31 

Value y (100) , when computed on one machine having 10 decimal preci- 
sion will have + — error bound. On the other machine having 

1 

20 decimal precision computable number y(100) will have + Jqq 
error bound. In Section 4, we will see that two computable number 
are called equal if their graphical representations have a common 
overlapping region • Here value y has to be in both the computed 
numbers. Thus we can say that both machines will give the same 
computable number. 

Difficulty (3) is that sometimes MAP real analysis 
results are totally unacceptable . In computable analysis, the 
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machine will give a result with the required error bound, 
although this may require many mors computations than MAP . 

If result cannot bs obtained to within required error bounds f 
in the programme provision can be made to indicate the situation. 
Thus there is no chance of unacceptable results going unnoticed. 

Difficulty (4) is lack of control on error in MAP 
real analysis cannot serve our purpos e because there is no way 
of computing result with the required error bound. Computable 
numbers by definition assure about that. 

3.2 A Sub-analysis of Real Anal ysis 

To understand the difference between MAP real analysis 
and computable analysis further, it is useful to ask the question 
1 Is computable analysis a sub— analysis of Real Analysis? 1 There 
are two ways to look at the relationship between real analysis 
and computable analysis. 

One way is to start from the fact that computable 
numbers are real numbers but n4s(^all real numbers are computable 
numbers. Yet, a good body of results of classical real analysis 

holds for computable number too. In this sense one can say that 

* 

the computable analysis is a sub— analysxs of real analysis* 
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Another way is to include in the number representa- 
tion the error bound t within which a number is known. This is 
the computable number concoct. Real numbers can be taken as a 
limiting case of compute ole numbers with error bound reduced 
to zero. In tnis sense one can say that real analysis as a 
subanalysis of computable analysis. Following Aberth 
we prefer the first approach. 

3.3 Basic . Interconnections 

Ejqjo^riag differences between MAP real analysis and 
computable analysis we can recall various interconnections of 
these both with error analysis and numerical analysis. An 
overall idea of these interconnections is pernaps best had 

r 1 

from the diagram shown in figure r .As sho'.m the re t 

% 

the MAP real analysis does not provide error information but 
only the value. Interval analysis is the name of the analysis 
which carries MAP real analysis and error bound computation 
both while computing. Computable analysis goes one step ahead 
by regulating maximum error to the required value by proper 
feedback procedure. 

3.4 Histori cal .Comments 

The full picture about computable analysis would not 
be complete unless we take into account certain concurrent 
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developments in constructive analysis and interval analysis. 

Cojistj^cl^ ; Constructive analysis can be 

said to have emerge^ in concrete form from the work of Bishop 

. In his book 'Foundation of Constructive Analysis' he 
carried out a reformulation of real analysis in terms of the 
concept of functions and numbers that can be computed. A more 
recent book on this analysis is by 3ridges. Stress in constru- 
ctive analysis is on what the results of important real analysis 
theorems would look like once we follow the constructive approach. 
The stress on computation is however, only to the extent that 
the notion of error bound is incorporated in the definition of 
a number. But the question whether a program can be obtained 
to determine numbers within given error bound is not dealt with 
at all. 


Despite its limitations, a positive contribution of 
the constructive approach is in bringing about the recognition 
that the phrase ‘it exists 1 better be interpreted as ‘it can be 
computed ‘ . 

3.4,2 Interval analy sis, as the name also suggests, takes 
intervals of real numbers as new kinds of numbers. In it a 
number is replaced as the basically by a pair of oral numbers. 

[a, b] . The numbers a‘ and b‘ are refered to as left and right 
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end points. This indicates tnat computable number is also a 
parallel concept to interval analysis, where 

[ a, b ] = [ £(s) - ~ , X(s) + i ] 

for the com, 3 a table number X(s) representing a real number 'X 1 that 
has been reorcsented as a>b in interval analysis. 

In the literature [6] we find two reasons of 
interval anthmatic : 

1* Sxa c__t _I n_t e rya 1_ Ar_ijhxaa tp c: s In which, we define operations on 
interval numbers as follows s 

[a/b] + [c, d] = [ a + c, b + d] 

[ a, b ] ~ [ c, d ] = [a - c, b - d] 

[ a, b ] . [ c, a] = [min (ac, ad, be, bd) , max (ac, ad, be, bd)] 

and [a, b ] / [ c, d] = [a, bj . [ ^ , ”] 

2* Rou nde d Interv al, A ri thmetics s This includes a practical way 
to find error bound, by taking differences between single pre- 
cision and double precision results in interval analysis. 

An important fact about interval analysis is that the 
Distibutive law does not hold for worth operations. 
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Let us see this with this example : 

[ 1 , 2 ] ( [ 1 , 2] _ [ 1 , 2 ]) ft [ 1 , 2 ][ 1 , 2 ] - [ 1 , 2 ][ 1 , 2 ] 

R.H.S. = C 1, 2 ] ( [1, 2 ]_ [1, 2] ) 

* [ i # 2 3 [ -l, i] 

"C"*2 / 2 ] 

L.H.S. =[ 1, 2J[1, 2] - [l , 2] [ 1, 2] 

--=[1, 4]-[l, 4 ] 

= [-3, 3 ] 

Thus we have relation 

[1, 2 ] ( [1, 2']- [1, 2 ])C[1, 2][1, 2] _ [1,2][1, 2 ] 

3.4.3 ^rnpujtahLe Analysis 

Aherth's book 'Computable Analysis' is the collection 

of work in this direction and this book is the mam source for 

this thesis also. In this book a computable number is defined 

with 

as a finite step program whict/ finite computation can produce 
the needed result with the required error-bound on it. Let us 
say we want to compute 

y = f <x) 
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Then computable analysis will give as an input, the 
error bound on y and tha computer program has to decide about 
the error bound on X(s) . With this changed eoncep t of error, 
computable analysis incorporates the constructive idea as it 
uses a program as a primitive element and computability on a 
finite stace machine as a proof for any compu'cabl 3 number. 

3.4.4 Interval, Analysis versus Computable Ana lys is s After this 
historical review, we can think about the relation between interval 
analysis and computable analysis as interval analysis has been 

tried out on computer [f] rather extensively. It would 

« 

seem the algorithms and software techniques developed for interval 
analysis can be usefully employed in computable analysis. 

In interval analysis, every number is represented by 
an interval say a, b for X such that 0 . < x < b. For calcula- 
tion of a function, independent variables are given as inuervols. 
Calculation will take care of these intervals and resulting 
interval number is computed. This resultant interval is not known 
in advance. 


Computable analysis on the other hand specifics the 
required maximum error bound on the result. In figure 3*1 (a) 
we have seen that computable analysis have interval analysis as 
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prj.rtj. so the work done in lncurval analysis area can help 
compute ol - analysis. One suen atuempt will be made now. 

3,4.5 Considerations “ J.I1. '/aka in his paper 

‘Portable software for ixiterval a rithmatic ’ [7] / has 

indicated that three different versions of Interval Arithmetic 
Software are running on computers. Ja will enumerate some 
features which can help m developing the software for computable 
Arithmetic 

(1) Brent's^ mul^tpole^jpre^cision^ pp.ckacje s In this package an 
initialization routine sets paraneters as base precision and 
exponent range. Then it allows a programmer to use multiple 
precision but intermediate precision should be less than maximum 
precision fixed while initialization. This is aporently to 
carry out a chain of computations on a variable precision than 
the present fixed mode. This facility can be used in compu- 
table number programs as error-bound of inpurs is varied till 

we get the required error bound on the result. 

(2) Single precision Interval Arithmatic package s This package 
was interfaced with Brents multiple precision package. 

This idea also can help in preparing corrputable ana- 
lysis software in modular form. First module can give programs 
for standard functions as per need of the computable analysis. 
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This module can be interfaced with some variable precision 
module, analogous to the Brent's multiple precision package. 

(3) Interval software with D irected r ounding s interval software 
with Directed rounding for interval arithmatic was implemented 
by writing software to perform floating point operations in 
machine language (preferably) . Although not much detail is 
available about Directed rounding is provided but it seems to 
be meant to provide a facility to change rounding errors in the 
program itself. 

(4) Nonstandard V ariabl e ; About precompiler they have indi- 
cated that Nonstandard variable were declared to specify interval 
number and precompiler is to generate standard Fortran programmes. 
In these operations and functions on nonstandard data elements 
are replaced by CALL'S to needed modules of supporting packages. 

This idea is also useful for the computaole software 
as nonstandard type elements here will be 'computable numbers*, 
in place of interval numbers. 

More specific comments cannot be made at present, but 
in future these softwares of Interval Arithmatic can be obtained 
and examined for the above mentioned facilities. Computable 
software can be prepared on analogous lines which will make the 
use of computable analysis practically feasible. 



Section 4 


Th e C omputable, . Numhe_r_ JField 


4.1 In cr od.ucti.on 

From the discussions of SectionS2 end 3/ it should 
he sufficiently clear why real numbers need to be replaced by 
computable numbers. We shall now, in this section examine the 
mathematical details about computable numbers. We shall in 
particular examine the fact that like the real numbers, com- 
putable numbers also constitute a field. We shall see how the 
various anthematic operations and relations for computable 
numbers are to be interpreted. In the last, concepcs oj. computa 
ble functions, computable sequences and complex computable 
number will be discussed. 


4 . 2 The .Basic ^Concept 

While talking about computable numbers in Section 1 
we had broadly explained it to be a programme X(s> that pro- 
duces fron the value of s given as input, the value of the 
number X to within an error of + 1/s. 

Now, what is the relationship between X(s 1 ) and X(s 2 ) 
the programme outputs for two different inputs s 1 and s 2 ? 
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'To explain this, let us take a programme u(s) that computes the 
value of 71, Let us say X(50) is to be computed, wmch means 
that the value of u is to be computed with maximum error + 1/50. 
This is shown in figure 


J.._ J 

1 i, 


= __„t 


- — - 


X(50)-— X(50) X(50+e~) : 3.11 


3.13 


"1 

3.15 


50' 

pig. 4.1(a) 

Likewise to coraoute X(100), the result should have 
a maximum error of -r 1/100. This can be represented graphically 
as in fig. 


X(l00)-^~r X (100) X(IOO)-}™ 


1 3.13 




— — j- 


3.14 


3.15 


Pig. 4.1 Co) 

’I 

Together, X(100) and X(5J) be can be represented as 
shown in Figure 


X(100) 


'loo 


X(10 °)-Too 


L 


X(50)-ji XCO) X(100) X(50)-l— 1- 

pig. 4.1 (c) 

Here we see that X(100) and X(50) have an overlapping 
region (dashed) • Let us see what are the possible ways in 
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which the overlapping can take place. There are essentially 
three different v/ays as shown m figure 4.2. 

Tor the number Y = 3/ and its confutable values 

Y(50) and Y(100) may be inucuelly disposed. 


Case 1 


Y ( 100 ) - 100 

JgJ — 


Case 2 


Y(50) — ~ Y (50) y Y(50)+ ^5 

F±g. 4 . 2 (a) 

1 


Y(100 ) ?TS5 
•( 100 ) 





Z(50)_5 C 


t (50) -h 


50 


|— 


Z (100) S(loo) 2(100 ) -^35 2 ,( 50 ) 

Fisr. 4.2(b) 


Case 3 


4 


X 1 (50) -t 


50 


X 1 * 5 °) ' i 50 


x 1 (loo) Xj. ( 100 ) Xi ( 100 ^Too x i4oo) 

pig. 4 . 2 (c) 

It can be seen that it is not possible for X^O) 

and X 0 (100) to be such that 

£* * <1 


X 2 (50) — X 2 (100) < 100 ■' 50 
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For, it were so then 

X 2 (50) -5*d X 2 (50) ' : 5§ 


X 2 (100 >"l|o X 2 (l00 > -2 (100) - , T00 X 2 (50) 

4-2(d) 

But as shown in figure 4 . 2 (d) the intervals of (100) and 
X2(5G) do not overlap which is not possible if X 2 (100) and 
Xo (50) are both to approximate the same number X. 

Witn these clarifications, we are now ready to 
look into the precise definition of a computable number. As 
wo presently see one first needs to introduce the notion of 
programmable function followed by that of a computable process 
and the equivalence of computable processes. Our presentation 
is based entirely on Aberth [ l] . 

Let P denote a programme on a URM and let P(x^,...3c^) 
denote the output of P when. P halts. If P does not halt then 
P(x 1 ,...x n ) is said to be undefined) , where . . ,x n are values 
of the programme input variables at the starting point. 

Let <P(x L ,..x n ) be a rational valued function of 
the rational variable x x ,..x n and let p bo the corresponding 
programme such than P(x jL ,..x n ) is defined if and only if 
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( P(x 1/ ..x n ) is defined, we say than C is a p ro cj r amma b l e fun c t - 
ipn if 

P (x 1 , . . .x n ) =¥ (x 1 , . . ,x n ) 

The programme P is said to reahize ■ . 

A C CTno utablc^ process is a programmable function a(s), 
that is defined for all po stive int-gers s and is such that 
for any two positive integers s and t 

j a(s) - cc(t)| < ~ ^ ( 4.)) 

wc may note at this point that the justification of the inequality 
(■4*1 ) lies in the arguments of the overlapping regions given 

at the beginning of this discussion. 

Let a and e 1 be two computable processes, we say, 
that u and u‘ are equivalent, u s. if for any positive 
integers s and t 

| u ( s) - U« (t) | < | | ( 4 * 2 ) 

Condition 4*2 states the equivalence relation of computable 
processes a and u l , which is reflexive transitive and commu- 
tative in nature. 
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We may now finally define a computable number as 

follows. 

The equivalence cl_iss of e computable process is 
a i^^Jitablp^ number. 

For the sake of notation convenience we snail use 
X(s) to define a computaola numoer approximating X. 

4.3 Error, Representation in Corrnutable, lJurrber 

So far we have interpreted X(s) to mean the value of 
real number X to within an error 1/s. Thera are other equally 
useful ways of describing error in X(s), For instance X(s) 
is a computable number then we may choose to define error 
equal to 4 log^s, Likewise on, instead of specifying s, may 
specify E = 1/s. All these choices are matter of notation 
convenience only. 

4 . 4 Th e ineq uality. Relation 

To be able to fully define the arithematic of 
computable number it is now essential to formulate various 
arithematic operations we first take up the inequality re- 


1 atio n, 
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Let us see some examples to understand the notion 
of inequality in computable numbers. 

Let X 2 (s) be a programme to compute V" 2 as given 

below : 

1. value of s given 

2. n = 1 

2 

3. (n/s) - 2 if positive, n/s is the result 

otherwise go to step 4. 

4. n - n+1 

5. Go to Step 3 

Let us say "^(s) is the computable number approxima- 
ting V 2 produced by the algorithm given below 

1. Given value of Sj 

2. x = o, r = o 

3. n = 1 

4. (X-i-n.lO r ) -2.0 if negative n - n-KL, GO TO 3 

1 ? 

otherwise x = x t(n-l) 10 
x - x - 1 

GO TO 3 till x = -s^ 

S 1 

(Here x - 10 ) * 
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Now let us search if X 2 (s) is greater than, equal 
to or less than Y 2 (t) for ‘s' and 't* positive integers, with 


tnese results 


X 2 (50) 


X 2 (50) « 1.42 

t 


Y 2 (100) = 1.41 

— ! 

1.40 

1.42 

X 2 (50) = 1.42 

Y n (1000) - 1.414 

I 

Y 2 (100) 

2 

1 

1.40 

1.41 


| 

y 2 (1000) 


1 

1.413 

1.414 


1.44 


1.42 


1.415 


As we see from the diagrams, since X 2 (s) and Y 2 (t) 

must contain the number X, that two will always overlap. When 

for any arbitrary values s and t, if such overlap exists for 

X (s) and X_(t) then we say they are equal otherwise they are 
A £ 

unequal. 

If there is a choice of some s and t such that over- 
lap does not take place, then we say two computable numbers 
unequal. In particular if X 2 (s) and Xj(t) are two computable 
numbers then we say XjU) < Xg(t) if for some choree of s 
and t,X 2 (s> a | < X 3 (t) - | . rfe can alternatrvely say 

that X^Ct) > X 2 (s). 
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4 * 5 Operation on .Computable lumbers 

Let us no v/ turn to arithematic operations on corrput— 
able numbers. First we need to define arithematic operation 
on computable processes. 


c</ and p are computable processes and operation 
l 0‘ represents one of the rational operations -h, i and % 
then o,Op is the process given by 

(0.0/3 Ms) = u(s- ) 0 p(s 1 ) 


where 


min v {max I t(v)op (v)-u(v) t ~)0(/3(v) ~)J 

v > 1 "" v ~ v 1 


s / 


It can be shown that uop is also a computable process. 
Here division operation holds only when P ( ^ 0. Thus opera- 
tions between two computable numbers are related with error 
com ju tat ion and control metnods also. 

Definition s Let Cu denote the computable iwe. the equiva- 
lent class of a computable process a. Let ‘O’ denote one of 
the operations +,-<* and then a p is the computable process 
c„0£ . For divide operation, to hold we must have £ (s^ A> 

for any positive integer s-j_. 
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On the same lines, we Gan reformulate the concepts 
for computable functions and computable sequences. 


4.6 Co m p u tab le ^functi on a nd p^ejquence 

Considering the computable number as a programme the 
concept of function has to be recast m cerms of programmes and 
computable processes. 


Consider as an illustration the following function 
in real analysis. 

2 

F(x) = X H- X 1-1(5 


Here x is represented as S (x) . the computable numbers appro- 
ximating X. further let Y{s) be the program to compute f 3 
with + 1/s error bound. 

Then the computable function corresponding to F(x) 
F(X(s 1 ) / s) = X(s 1 ) -i-X^Cs^ + 


wh: 


ire 

1 

s 


max 


! X(s 1 )- r X 2 (s 1 )+Y(s ;L )+ | )-(X( Sl )i | ) 2 -Y( Sl )± | _>| 


or 


Here 


1 2 , 2X(S 1> , 1 

PP 7— J 

S S-j^ Si S-j_ 

s = min v{mox|F(x(v) ) - F(x(v)h_ - )|- g'* 
1 v 1 



Thus error control method is just parallel to used for the 
operations on computable numbers. Sequence when used with 
computable numbers will have changed formulation. This we 
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will see now. 

In engineering or scientific work sequences occur 
very commonly. A precise definition of computable sequence 
is given as follows, 

Let e n (s) a programmable function, when n is a 
positive integer, such that for a fixed n, cc n (s) is a computable 
process. We can say that u- n (s) defines a sequence a n where 

~ u n * 

It can be shown that if a n is a sequence that if 
is a sequence then so are the and given by 

n m 

b -- S a. an<3 - c -n ~ n a i 

n x n ±=1 x 

This result is of basic important in the analysis 
of digital filter algorithms. 

4 . 7 Comp lex Co mgu_teble^^al^ir 

Once we replace real numbers by computable numbers 
we can extend the notion of .computable' in complex field also 
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complex numbers can be replaced by 1 computable complex numbers ‘ » 
And then parallel to computable analysis, a computer cornpu cable 
analysis can be x/orked out. 


A complex number is given by a pair of rsal numbers 
*' a l /a 2^ w ^ ere a j_ is called real part and a. 2 ~ as imaginary part, 
A computable complex number is given by a pair of comoutable 
numbers (A-^^) as in a complex number, here A^ is a programme 
^1 ( s l) ' as r3a -^ Pact and A 2 - another programme e^Cs^) as 
imaginary part. 


Now we come to most crusial result namely that 
the computable numbers together with operations defined above 
form a field. This is known as Rice Theorem [ 1 ] . In viex/ of 
its importance, we x>r3cent here a structural version of the 
proof of this theorem. 


4 . 8 Structured Jiroof _of Rice- -Theorem 

The structuring of the proof is based on the guide 
lines given by Leren [ 9 ] to make it mors readable. The 
basic idea of structuring is as follows, proof is broked into 
different levels and each level will be sole to convince the 
reader, what ne^t is needed to prove the theorem. Vie will see 
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this with Rice-theorem' s proof 

Leycl^l - Theorem is computable numbers form a field. 

bev el 2. ~ In level 1 1 to prove that comoutaole numbers form 

a field/ we snould try to prove that computable 
numbci.s fullfil all conditions necessary to call 
it a field/, these conditions are that computable 
numbers obey 

Pi Commutative Lav/ 

P2 Associative Law 

P3 Distributive Law 

P4 identity Law and 

P5 Inverses exist 


Go now in next levels we need to prove that above 
coriui tints are satisfied for the computable numbers. 

Lov'd, 3 c To prov^ it to P5 propositions or conditions we will 
try to tee for ouch proposition/ how to prove it and what 
tluns will bo need to prove it. 

For (pi l j To prove that computable number obey commutative law 
we have to prove that computabl e numbers obey these relations 2 

Pl.l ally * b + a 
PI. 2 a. x b, *= b x a 

wher 2 a,& b. both am conputable numbers. 
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Pojc ^(P_2_)^ s If associative law holds good for the computable 
number or not, can be known if we can pro/e that comsutaole 
numbers obey these two relations; 

P2.1 a-f-(bvc) = (a-rb)-fc 

p2 , 2 a x (b x c) = (a x b) x c 

where a, b, and c are computable numbers. 


F or P3 z For computable numbers, if one can prove that ror 
two computable number a and b 
a(b-fc) = ab + a.c. 

holds then distributive property for computable numbers is 
satisfied. 


Prom_P4 s 
el aments. 


For proving that cortnu cable numbers have identity 
rolations 


P4.1 A-i-0 = a 
P4.2 bxl = b 

should hoi a where a and h are computable numbers 


Fro P5 : For the last property one has to prove 
table numbers have their inverses, m other word 


that coirtpu- 
relations 


a + b = 0 
a x b = 1 


for a and b computable numbers 
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iiSZS.'L.-A ° -*- n level 3 we Game to know, what is needed to 

prove the proposition of level 2, Now we will see each propo- 
sition and will prove conditions required (hy level 3) for 
computable numbers. 

Fo r P i s From level 3 relations to be proved are 
P 1 « 1 a bb — b *i* a 

Pi. 2 axb=bxa 

where a and b both are computable numbers. 

Pl.l 

Relation is 

a -r b = b a. 

As a and b are comou table numbers, let us represent these by 
a(s) and p (s) computable programmes, so 

(a + J3) (s) = t( Sl ) + 6(s 1 ) 

(p -I- a) (s) = p (s 2 ) '• «(s 2 ) 

To prove (t -!- p) (s) = (P -l- u) (s) 

we need to prove 

u(s^) + P (s^) — P ( s 2 ) + t ( s 2 ) 

We choose Y value such that it is m overlap region 
of ^( Sl ) and a(s 2 ) and «Z» value from, overlap region of p (s^ 


(4-1) 

( 4 - 2 ) 


I 
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and 3(s 2 ) and putting these values in ( 4- D and (4.2) equa- 
tions . 

^(s 1 ) + p(s 1 ) = Y -i Z 
3 (S 2 ) -i- a( Sl ) = z + y 

Since Y + Z = Z + Y for normal plus operation. 


P1.J2 : On same lines 

(a * p) (s) = a(s 2 ) *' p(s 3 ) 

(3 /v o<) (s) = B(s 4 ) * o.(s 4 ) 

Let M and N be two values from overlap region of a (s^) / o,(s^)' 
and p(s 3 / p (s^) respectively 

a(s 3 ) * 3(s 3 ) = M * N = N * M = p(s 4 ) c.(s 4 ) 


ffftrP2 t 
relations 


To prove Associative law for computable numbers rea 

to be proved are 

p2«l a +(b+c) = (atb)-hc 

p2 ,2 a x (b x c) = (a x b) x c 


where a^b/C are computable number^ which can put as cc(s) / J3 (s) 
and r(s) respectively. 


P2.1 Relation is 

a -i-(brc) = (a+b)-i-c 
(a-i-(3+r) ) (s) = ((u+p)+r)(s) 


or 
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Using (Ou.0/3 ) (s) — ^(n )Op(u ) relation. 

s s 

RHS (tc+(£H-r) ) (s) a a (s 1 )+(p+r) (s 1 ) 

= t(s 1 ) b(s 2 ) + r(s 2 ) ( 4.3 ) 

LHS ((u+p) + r) (s) = (Cv-;-j3) (s 3 +r(s 3 ) ( 4.4 ) 

= Cv(s 4 ) + $ (s 4 ) v r(s 3 ) 

Taking P value from overlap region of a(s 4 ) and 0 .( 5 ^), 

Q value from overlap region of p( s^) and p (s 2 ) 

and R value from overlap region of r(s 3 ) and r(s 2 ) 

From equations ( 4*3 ) and ( 414 ) 

a( Sl ) + e(s 2 ) 1 - r(s 2 ) = P Q R 

<x(s 4 ) + 0(s 4 ) i- r(s 3 ) = P -r Q -i- R 

Hence (a+< p -!-r ) ) ( s ) a ( ( 01 -p ) -i-r ) ( s ) 

P2« 2 _s ax(bxc) = (axb)xc 

or (u*p ,f r) ) (s) = ( (u-' f p)“r) (s) 

On the same lines . 

(a* (j3*r) ) (s) .» t(s 1 ) Vf (P"' : r) (s 1 ) 

« a(s 1 )* 3 (s 2 )-r(s 2 ) ( 4 * 5 ) 


((a*0)*r(s) = (a*p) (s 3 )*r(s 3 ) 

= tt(s 4 )* 3 (s 4 )*r(s 3 ) 


(4.6) 


! 
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Selecting P # Q,R as above mentioned - 
From equations 4- .5 and 4 . 6 
u(s 1 ) A j3 (s 2 ) <f r(s 2 ) = P'-Q'-R 
u(s 4 )*p(s 4 )*r(s 3 ) - P V 'Q 'R 

Hence (c.* (p+r) ) (s) = ( (<-*3)* f r) (s) 

For P3 : To prove Distributive law for compucable numbers; 
relation to be proved are s 

P3.1 a (b+c) = a.b + a.c 

where a/ b and c are computable numbers representable a(s), 
P(s)/ respectively. 

Relation to be proved is 
a (b+c) = a.b a.c 
or (u. (p-i-r) ) (s) = (u.p -t '■-.r) (s) 

(&*(p+r))(s) - c^Cs-l) * (P’rr)(s 1 ) 

= u(s 1 )*0(s 2 )+r(s 2 )) 

= a (s 1 )*P (s 2 ) + t(s 1 ) w r (s 2 ) 

(&.*p+o.*r) (s) = (a*p) (s 3 )t(c^r) (s 3 ) 

= u(s 4 ) A 'P(S 4 )H-t(s 5 )“F(s 5 ) 


(4.7) 


(4.8) 
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Selecting P value from overlap region of u.(s.j_)u(s 4 ) and a(s^) 

Q value from overlap region of £ (s 2 ) :nd p (s^) 

and R value from overlap region of r(s 0 ) and r(s^) 

From equation (4*7) and (4.8) 

u ( s 1 ) "'p ( s 2 ) +a ( s-j^ ) -'r ( s 2 ) = P 'Q-.-P-'R 
a 1 s^) ,f p (s^J-l-u-Cs^) c r(Sg) = P-'Q+P R 

Hence (a* (j3+r) ) (s) - (u*p K- ,e r) (s) 

For P4 : To see if identity elements exist for the comput- 

able number set# relation to be proved nee 
P4.1 - a -i- 0 = a 
P4.2- a x 1 - a 

where ‘a 1 is a computable# which we can represent b_ t (s) :# 

Here (a 1 0) - ^(s^) ~' r 0 
= u(s-j_) 

and (a*l) (s) = o-(s 2 ) 4 
= u(s 2 ) 

For P5 s Property of 1 esistence of inverses 1 can be proved 
by the relations 

a b = 0 

a x b = 1 



where a and fc> are two comp unable numbers and we represent them 
as u,(s) and |3(s) respectively. 

P5-.1 i a + b - 0 

or (u+p) (s) ~ 0 

or o-Cs^) + p(s 1 ) = 0 (4.9V 

P5J2 : 

a x b = 1 

or (o.*j3) (s) = 1 

or a(s 2 ) ; 'j3 (s 2 ) = 1 (4.10) 

Seeing results (4.9) and (4-10) , we can say that compu- 
table number set fulfils property of existence of inverses. 


§1867 



Section 5 


How JTo Make It _ JP ra c ta cal 


5*1 Introduction 

To bo abl o to make use of computable analysis for 
computational purposes, we have to look for practical approaches 
also. But it must not increase the complexity of programme for 
solving numerical problems at the same time computational algo- 
rithms based on it should automatically provide an indication 
of the error bound at eacn stage. 

Some of the modifications n?eded for the purposes 
of general computable number programs would be the following 
ones s 

1* Library functions will have to be rewritten. 

to 

2. operations +, * and —will have/be programmed In 

computable number sense. 

3. Programmes for the comparision of two computable 
numbers as well as for binary relations in general, 
will have to be reformulated. 

. Infinite series expansions as trigonometric functions 
will have to be constructively redefined. 


4 
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5. In working backwards from the accuracy of the result 
of an algorithm operation/ an appropriate rule will 
have to be evolved, for distributing it to operands. 

6 . Modularity in programmes will become all the more 
important. 

7. In view of the increased computational work that is 
required to incorporate error as a parameter/ method of 
cutting down the number of iterations required for 
computing a function to within a given error bound 
will have to be devised. 

Let us now discuss the need for each modification and some po- 
ssible ways to carryout that modification. 

5.2.1 Lihraix Actions 

Library functions on computers allow us to compute 
standard functions with most accurate possible results at 
present. For computable analysis we have to rewrite these 
library functions^ Such that the programme should be able to 
compute the result for the given final error bound. 

For this, solution lies in arranging for the compu- 
tation of error bound and error controlling feedback. Let us 
illustrate this with the help of an example. Let us compute 
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0.5 


with + 


100 


error hound. 


0.5 


Now e Gan be put in its expanded form as follows 


e °' 5 = 1 + 0.5 + + I°-ii + ... + 10,52 

i2_ 13- UL 


n 


+ . 


For having + error bound, tnis series has to be 


th 


truncated at the n term such that 


100 


< - 


(0.5) 


n+1 


m+1 


IP.'AL 

sn+2 


n+2 


+ ... 


Thus algorithmic stops needed are as follows s 

S "t 

1 . start summation from 1 term 

2. compute next term 

3. Add it to the sum 

4 . Find error due to truncation on that term, 

1 

5 . check if truncation error 

6. if not go to step 2 

el.se stop execution. 

In this procedure truncation error has to be computed. 
One way to do this is as follows 


E 


(0,5) 


n+1 


t i . ° « 5 ... wftZ -i- 't 

( 1 r n+2 h (n+2) (n+3l '*• } 


(0.5) 


n 


(n+1 


i0j5± 
in hi 

T i > . n i 


ri+l 


(1 + r. + * r 2 + r i * r 2 * r 3 + • • * ) 
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where 



0.5 

n+3 1 r 3 


a*5 
n+4 * 


* * « ^2 t C « 


As we have ... for n, a positive integer, by putting 

r l ~ r 2 ~ r 3 WG can wits 


n — 


(O.t) 

*JT- 

in-rl 


n*l 


( 1 + r. + 


or 



Xp.5) n+1 

Tn-t-'f 


1 



(as < 1 for a converging series) 

( 0 r) ^ 

This -h value can be taken as the truncation error 

| n-fi l-r. 

j. "tJhi 

by truncation of series at n term. 


Such algorithms have to be written for each library 
function. Once such library functions are available, computable 
number programmes can use them directly by calling these 
library functions when ever required. 


5 . 2*2 Operations, 

In section 4, we pointed out that for the computable 
numbers an arithematic operation is Of the form 

Z(s) = X(s 1 ) 0 Y(s 1 ) 
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where X(s^) * Y(s^) and Z(s) are corrputable numbers and *0 1 
denotes an arithematic operation, ^urthar* 8^ is related to 
s as 

Si = min v { rrraxl X(v)0 Y(v) - (X(v)-1)0 (Y(v)-i)| <- } 
V* > 1 “"V *“"V s 


This means that an arithematic operation in compu- 
tet 1 analysis has to be treated differently and not as we 
carryout an arithematic operation on two MAP real numbers. 

To bo more specific, we will have to include in the program 
implementing that operation a provision for obtaining from the 
value of s the error in the final, the error to be allowed in 
the operand of the intermediate steps. To give an example, 
consider the arithematic expression 



x 

3 


In terms of computable numbers, this becomes 

5 *x 2 (s 2 > x(s 3 ) 

Y(s) = xtsp + j -3 + — 3 

Suppose it is given that s = 100 and X(s^) gives the value of 

real number Sin 40° with + — error. To compute Y(s) we can 

1 

proceed as follows • 
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Here, there are two binary operations involving three 

terms. Let us distribute the error — equally. This 

s 

means that every term should be computed with + ~^ ’ q 
error bound. 


2. Taking each term one by one. 
(a) X(s^) will be computed for 


1 = _ 1 
300 


5 . 2 

(h) * X (S 2 ) will be computed using computable 


number algorithm for the operation • x, # choosing 

( =! ) is 

13 A 1S 300 


So such that total error on ^ * X (s 0 ) is 


1 /2 

(c) Swuaro root (XCs^)) is computed using library 

function programme ‘square root 1 such that s-, bb 

1 X( ®3 >-5 

selootoa to havo the error bound 355 on 3 — . 


This way computable operations will have to be carried out 
using a separate algorithm which can control error to the 
required value. 


5.2.3 Comparision of two computab le n umbers 

For cortparmg pwo computable numbers, we have to 
call some subroutine to decide whether one computable number is 
greater than or equal to or less than the other computable number. 

Let us, say X, Y are two real numbers represented 
by computable numbers X(s) and Y(s) respectively. We c a n have 
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three possible cases i 

1 1 

1, Xf X(s) + ~ < Y(s) - ~ , then X(s) < Y(s)^ this can be 
graphically shown as in figure 5.1(a) 


p &*- MMIli * -1 T . ||»>. hm L , J, 

! X(s) ! 


-r i 


X(s). 


X(s)-ti Y(s)-» 

o S 


Fig. 5.1(a) 


Y(s) 


1 


Y(s) 


2. If X(s) ~ ~ > Y(s) + ~ then X(s) >Y(s), this can be 
s s 

graphically represented as in figure 5.1(h) 


1 Y < s ) ^1 \ x(s) ^ 1 

Y(s)~~ Y(s)+^ X ( s ) X(s)+^ 

F%. 5.1(h) 

3. If X(s) +| > Y(s) - ~ and X(s) - ~ < Y(s) + - 

then X(s) = Y(s) 

Tnis is shown graphically represented as in figure 5 ,J> 1 ( c ) 




4— ••"y(¥)' 

X(s)t| 



Fig. 5.1(c) 
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Thus all three inequality relations have to be tackled in 
the new sense. 


5.2.4 Trigo no metric fu nctj-ons 

Trigonometric functions on general have infinite 
terms in the expanded form. We will discuss here a case of 
convergent series with infinite terms and see what the problem 
is in computing error bounds and we will try to see if there 
are reasonable solutions. 

To understand the problem, let us consider an example 

to commute log(l+X(s), where a real number X xs represented by 

X(s), a computable number. An expanded form can be as follows 

2 3 

log ( 1 + X (s) ) = X(s) - — ... 

In such a case there are two variables in error 

1 

commutation. The first one is error bound + — , on X(s.) as 
per the computable number concept discussed in Section 1, 
Another one is the number of terms ‘n* after which truncation 
is done. So the problem is to decide the trade off between 
these two errors such that the summation of both is equal to 
the error required. In this respect first of all we will see 
tile contributions of each error separately for the function 


SinX as an example 
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For the function SinX, when we represent the real 

number ‘X* by a computable number X(s), the error due to 

1 

truncation and the error contributed by error bound + ~ or. 
X(s) hav 3 to be computed we will examine the error from one 
source with the other error assumed to be zero. 


Series expansion of SinX is as follows. 


3 5 7 

„ X X X 

binX = X ~ -rs" -i — rs- - 


X 


2n-l 


J3_ j5_~ tl ‘ *•' t2n-_l. 


. « • 


Co JL* whan « 0 or X(s) = X 
In the series expansion 


T„ 


n 


x 2n-l 

t2n-.l 


r 2n-l-l 


T 


n+1 | 2n+l 


T 


n+r 

T 

n 


X 

'2^C2nTlT' 


= r„ 


For sumrn t tion of first n terms. 


SINX = X 


3 5 ' 2n-l 

X , x , 

+ _. - ... + (2n-l 

1 — , U 

,2n+l 


X 2n+1 

i2nfl 


SM 


2n-l-l x X 

+ “j2n+r‘ 1 “ *f2n+2) (2n+l) + ‘(2n+2) (2n+3) (2n+4) (2n+5) 


+ ••• ) 


,2n+l 


= SM + ~]2SS. < 1 ■* r l + r l + r l + •" ) 
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T 


where 


n+1 


n 


T 


n+2_ 

t" 


n 


^2 

: 1 


since 


S r 


T T T 1 

n+l > n+2 > n+3 . , 

if “Tfr — - etc. so upper round of error 

n n+1 n+2 


X 2n+l 2 . 3 


! 2n+l 


( 1 + r^ + + r^ + ... ) = 


X 


2 n+1 


[2n-H * 1-^ 


(5.1) 


Here < 1 for series is convergent one. 


Case 2 » 


When the real number 'X* is represented by X(s) 

1 


with a finite error bound. + 


For this truncation error is 


taken to be zero, to see the effect of error bound + alone _ 
We have here the series expansion 

SinX (s) = X ( s) ~ -j- X-.C S A . . . 

(X(s) + |) 3 (X(s)jA) 5 

Sin (X(s) + -) - x(s> ± jj~ + —jg---- ... 


Thus 


E = max Sin (X(s)+~) - Sin X(s) 
s ^5 


l . 3xll§l„^X(sI+_l 


s 


2.s 


2s 


6s 


5X 4 (s) . lOX 3 (s) . 

' — I- J t 4. ( „ , 

120. s 120. s 

.. (5.2) 


represents the truncation error end represents 
the error due to ~ error in X(s). The total error E R can be 
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taken, as a first approximation, to be 


Equations 5*1 and 5*2 show that as the number of 

terms, n, is increased whereas E increases. On the otherhand 

s 

by increasing ‘ s‘ in X(s), one can reduce E and so that a 

s 

higher E T by summing fewer terms of the series may be 

This makes us think of many possible solutions for E and E , 

X s 

such that is constant. So the question of sel action of one 
pair of E^, and E g out of many has to be solved for minimum. 

Algorithms to compute infinite series that take 

into account this tradeoff between E m and E comoutable number 

T s 

programmes have to be devised, 

5*2,5 E qu al distributio n of error 

In operations with computable numbers of a distri- 
bution of error bounds to constituent operands is necessary, 

Z(s) = X(s 1 ) 0 Y(s 2 ) 

where '0* represents any of the operation +, * or-i,on 

computable numbers X(s^) and Y(S 2 )» Here the error ~ on Z(s) 
will be decided by and — ^ the error bounds on. X(s^) and 


l 
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Y(s 2 ) respectively. Several different pairs of s 1 and s 2 will 
produce the same value of s* Thus a well defined method has to 
ho prescribed to select one of many choices* Wo will not give 
an example to help evolve a method to select an appropriate 
s-j t s 2 pair* 


Let us compute Y(lOO) = V 2 + V 3 . This is particular 


of the 

form 




Z (s) = x 2 ^ s l^ 

+ X 3 (s 2 ). We will take 


^ 2 

= 1.414 ..., V"3 = 1.732 ... and V 2 -i- 

V 3 = 3.156 ... 

So we 

write some computable numbers as follows 


X 2 (50) =1.3 

X 2 (100) = 1.41 X 2 (200) = 1.415 

X 2 (SOO) = 1.412, 

X2 (1000) = 1,414 


X 3 (50) = 1.6 

x 3 (ioo) =1.73 

X 3 (200) = 1.735 

X 3 (0 

00) = 1.734 

X 3 (1000) = 1.732 


A few 

cases are now computed as follows 


(1) 

X 2 (100) +X 3 (50O) 

= 1.41 + 1.734 = 3.144 

(s=83, 3) 

(2) 

X 2 (200) +X 3 (200) 

= 1.415 + 1.735 = 3.150 

(s=100) 

(3) 

X 2 (500> +X 3 (100) 

* 1.412 -r 1.73 = 3.142 

(s-83.3) 

(4) 

x 2 (5o) +x 3 (iooo) 

= 1.3 + 1.732 = 3.032 

(s=47.6) 

(5) 

X2 (1000) + X 3 (50) 

=5 1*414 + 1-6 — 3-014 

(s=47.6) 
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So these cases show that lows* requires high s~ for 
the same s and conversely & low S 2 requiras a high A 

appropriate selection can he s-^ 55 555 200 

In general, for an addition of the form 

z(s) « x 2 ( S;l ) + x 3 (s 2 ) (5.3) 

in which the errors are related as 



it is appropriate to choose s^= s 2 - ^ . This choice is based 
on the argument that the error distribute should bring down the 
computations for each of the operands. For any number X(s), the 
smaller eh error " , the more the computations required. What 

S 

we want that for both the operationals X,(s^) and X(s 2 ) the 

11 1 

Si and s„ be so shared that while ~ = - 1- is constant, 

X 6 s 

Sjl + s 2 be minimised so that the computation for ( s.^ ) and 
X^(s 2 ) are minimised. This inplies that s 1 = s 2 « 

Consider now the multiplication m place of addition in eqn* 5*3 
Assuming X 2 ( s-j ) & ^(s£) and , s 2 large integers, 
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Thus here too ~ will require the 1 east compu- 
tational work. 

Here one question may arise* if this convention will 
held, good in general the answer is that except in seme very 
exceptional cases* it will always had good. Two exceptional 
cases are enumerated here, 

1, when X 2 (s % ) X^Csj) then s will be controlled by the 

larger value. 

1 X 3 (s 2> ^ 1 

s S 1 s i s 2 

2. Certain cases where X 2 (s) is a shorter programme thah 
X^(s) or in other word X^ (s) will require more computation nl 
work than that for X 2 (s). 

In such cases* one can choose the distribution of errors 

such that 

X 2 (s 1 ) x (s 2 ) 

«1 S 2 

5.2.6 Modula rity of Co mputable nsffitoe^PrgaaicTO 

Modularity of a programme helps in writing and also 
in understanding programmes. How this feature is useful in 
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the computable number programmes will be examined with help 
of three different cases. 

1. Let us conpute V2. This can be represents as a computable 
number X 2 (s) . Here the computable number programme will 
use * 2 1 and ‘c* as the input variables and will dive a 

l 

value after execution with + “• error bound. 

**** p 

2. Let us compute V2 $t 3 which can be represented as 

Z(x) = X 2 (s 1 ) * X 3 (s 1 ) 

Hero X 0 (s- ) and X,(s-) are two computable numbers of the 
typo mentioned in case (1). These ere used to compute a 
new computable number Z(s). 

3. Let us say we want to compute 

Y - Sin ( V*2 ' t (1 + ife ))/ which can be represented in 
computable analysis as follows 

y(s) - ^in (X 2 (s ± ) * ^ + X 3^ S 1^ 

Here X z ( Sl ) * (1 t XjCsp) is a computable number of type 

mention in case (2). That uses x^s,) computable numbers of 

type given in case (1). 

Z l (s 2 ) =X 2 (S 1 ) * (1 * X 3 (S 1^ 

Y(s) = Sin X(s 2 ) 


then 
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This Hay for given »s» s 2 should be decided but s 2 
itself is dependent on s 1# This case is most general 

° n ° * * In terms of which any arbitrary computation 
can be organized. More precisely any arbitrary corrpu cation 
will simply consist of a nesting of computation of the hind 
in several levels. 


Now we can see computational work in each case. 

C as e 1 

One can represent this as follows 

X 2 (s) a a(2,s) = SQUARE ROOT OF 2 

Here once ‘s' is specified programme X 2 (s) can be 
executed to compute the value. 


Casj3_j2 


This can be represented as 
Z(s) = a(x 2 (s 1 ) / x 3 (s 1 s) 

Hare *s‘ will be specified and programmes of X^s.^) and X (s^ 
progrommes are available. Let us see for s ~ 250, 

Z (250) = X 2 (s 1 ) * X 3 (s 1 ) 

where = min v { max | cc.^ (v) * Cu 2 (v) -(ct^ (v)+^p) ,v (a 2 (v)j ~ 250 ^ 
v> 1 


For v = 1,2,3, , . . 
-jL and put e x = 


we stop at when 

v. 


the error becomes less than. 
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Case 3 


This can be represented as follows 
Y(s) = Z 1 (s 1 ) * z 2 ^ s i> 


where s 1 -* min v I max 
1 v >1 


Z i (v) * Z 2 (v) - (Z 1 (v) i i)*(Z 2 (v) 2 i) 


< - l 
~ s 


For a given value of s start with v = l / 2,3 / 4 / ... and check 
the above condition on s 1 to be satisfied. That value of v 
gives s-j^* Hero r l\ (v) for each volue of v will ho 

Z x (v) = 2 ^( 82 ) * X 2 (s 2 ) 

1 1 | 1 

where s 0 ** min r [max |x,(r) * V r) - (x-r)^)” (X 2 (r)+j:)l \ 

r > 1 

Again X^r) for a given ar may be given as follows 
Xl (3?) = X u (S 5 )*3«1 2 (S 5 ) 

where s 5 « rain u { max| X^ (u)*X 12 (u)-(X 11 (u) +~>* (X 12 (u)4^| - ~ 1 
5 u > 1 

Evidently there is more computational work in case-3 
than in case-2 or case- 1 . 

This involved algorithmic approach can be made staple 
once we use modular programs for case-1. Programmes of case-2 
will use case-1 programmes. Programmes of case-3 type wall 
_ 1 and case-2 programmes. This modularity can be 


use case 




y&g # $*t now Ohavt off ® 
Mwab@r frograww# 

<2«e« l f(t) « V v) 

Oaw» 1 Fit) ® #1 it) 

« <^(0)0 


OaM 3 F(t) 01 


l 

I 



I 
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possible once we use the particular structure emerging from 
the flow chart of each case-1, case-2 and case-3, Using case-3 
one can cover the full, range of cornpu table number pro grammes, 
This fan ;r,\L si; rue Iran is shown in figure 5-2. 

5,2,7 Rejducjsioja in jZoj^dj^^ j fork 

The next question is whether we can reduce computatio- 
nal work. Let us see what compuratiouaal work; is involved and 
then wc will try to see if we can reduce it. Say we want to 
compute a computable number 

Z(s) = X 2 (s 1 )*X 3 (s 1 ) 

whore s x = min v ^ mQ;jc | x 2 (v) y ut 3 (v) - (X 2 (v) Ht~)j - } 

The approach to implement it, for a given value of ‘s' may be 

as follows 

1. start with v = 1 

2. compute X 2 (v) and X 3 (v) 

3. compute max |x 2 (v) ' V X 3 (v) ~(X 2 (v) " ( x 3 (v)i^) j 

1 

4. if value from step 3 is less than or equal to ~~ t 
stop execution else v ~ v+1 and jump to step 2. 

Here v is incremented by 1 in every iteration. Thinking about 


step (2) of the algorithm, if 
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X 2 (v) = Y 1 (u) *Y 2 (u) 

where u = mir^ v (max ) (r)*Y 2 (r)-(Y 1 (r) j~) * (Y 2 (r)^)} - i 

Thus for each value of v, ^ (v) has -to be computed and incre- 
menting v by 1 « We have to reach point when s-^ = v. The 
computational work is thus proportional to the value of s- L . 

To reduce the number of feedback loopings larger 
increments in v need to be considered, so an appropriate way 
can be to select a new value of v before going back for a new 
iteration, with the help of previously computed value in 
intermediate iterations. 


Let us say 

Z(s) * x 2 (s 1>* X 3 <s 1> 

where s, - min v { max } X 0 (v)*X, (v)-(X 9 (v)3-;) x (X n(v)t~) { <“5 


v >1 


Here max X 2 (v)*X 3 (v)~(X 2 (v) (X 3 (v)+~) 


X 2 (v) X 3 (v) 1 

— ' — — + i -*r 

V V V 


So for new value of v n we can use old value v^ and X 2 (v q ) 


and X 0 (v ) as follows 
3 o 


i Ml + v±. + 4 


V 


or 


«2< V o> + X 3 (v o> + -r > * 3 
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Here v n denotes value of v for the next iteration 
in feedback. Proceeding in. this manner can generally help 
cut down the number of feedbacks. We give below two examples 
of computations carried out in this manner. 

Firt one is regarding computation of 

z(s) ■= ^ r<t.i, si )*rci,2, gi ) 

whc>rc n = 1,2,3 

Second one is about the matrix inversion with com- 
putable analysis. Programmes used for both are listed in 
Appendix. In the result listed Tables (1) and (2), vre see 
that there is a great reduction in number of iterations on 
using suggested computation reduction approach. 
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Matrix 

A 



Matrix Inverse A 


3 1 

2 


0.3889 

0.0555 

-0.2778 

2 3 

1 


-0.2778 

0.3889 

O.Q555 

1 2 

3 


0.0555 

- 0.2778 

0.3839 

Value of i 

Maximum e, 

1 1 

error in elements of Matrix A, ~ = *-, 7 - 

00 7 

rror in elements of Matrix Inverse A/ 1/s ~ 

1/100 

Number 

of 

iterations (using Sl 

-= s il) = 367 


Number 

of 

iterations (using 

s x = X(s^s) = 2 


3 1 

2 

8 

0.3222 

-0.0111 -0.3044 

0.0400 

2 3 

1 

7 

—0.3444 

0.3222 0.0288 

0.0400 

1 2 

3 

9 

-0.0777 

-0.4111 0.3355 

0.0800 

9 8 

7 

6 

0.0666 

0.0666 0.0266 

-0.0400 

Value 

of 

error in 

. . 1 1 

elements of Matrix a, - * ~Y2sf 

1 1 

1 


Maximum error in elements of matrix inverse A, - = rdo ( SP-ven) 
Number of iterations (using s 1 ~ s^+1) = 1267 
Number of iterations (using s-j. = x(s^)*s) - 2 

Hore X(s n ) is the value of largest element of corresponding matrix 

JL 



Section 6 


Digital Signal Processing and 
Computable Numbers 


6,1 In tr oduc t ion 

Computing machines are used in digital signal proce- 
ssing (DSP), both for the design of DSP systems e.g. digital 
filters and for the implementation of such systems*. 

The software approach is now used extensively. Nume- 
rical algorithms and programmes have now beaome basic tools used 
in digtc'.l signal processing. The existing methods of realiza- 
tion and design of DSP systems are all based on real analysis 
in the maximum accurate possible (MAP) form. As a result 
they pooe difficulties that were discussed in section 3. On 
account of those difficulties, we felt that computable analysis 
would offer a more appropriate alternative as the computable 
number model specifically tabes into account the problems of 
machine representation. Thus we can think of using computable 
numbers instead of real numbers in digital signal processing. 

More specifically. 

. System theoretic properties of DSP algorithms can be 
reformulated in terms of computable analysis. 


1 
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2. Errors due to finite word length in a DSP algorithm can be 
looked upon in terms of computable analysis. In this manner 
the ‘Dead band effect' and the sensitivity of DSP algorithms 
can perhaps be better understood. 

3. Using the computable number concept can also be incorporated 
in the design of a DSP algorithms e.g. digital filters. 

The third task can be carried out in terms of computable analysis 
only when the practical tools suggested in section 5 have been 
fully developed. So we will discuss only (1) and (2). 


6 . 2 Ba^sic^rooertiej. 

First of all we shall look at some basic properties of 
digital filters to see how they can be rephrased. The proper- 

ti 2 S to be discussed are J 
1, Time invariance 
'2, Linearity 

3. Causality 

4, Stability 


Basically a digital filter Is a numerical algorithm 

. v,~ mrried out using available 

have finite confutations to he 

o -j-hr* already conrpufcod values of 
values of the input sequence ana the alruaay 
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the output. Each arithematic operation va wish to consider in 
the computable number sense. Thus we find a digital filter 
algorithm is essentially a computable process; 1 cd (section— 4) 
which produces output y(s) from input x(s 1 ) 

y(s) = a (x(s 1 )) 

whore 

s 1 « min v {max | a (x(v)) - a (x(v) + ~)| 1- } 

v £.1 -vs 

Here s in y(s) is the error parameter. The x(s) and y(s) have 
the computable number x(n,s^) and y(n, s) respectively as its 
value at the nth point. 

Time Invariance : Digital signal processing algorithm is said 
to be time invariant if 

y(s) = a (x(s 1 ) ) 
implies y T (s) = a (x T (s)) 

where x C (s^) Is the shifted sequence whose values are 
x T (a / s 1 ) = xCn-T/S.^) 

T . 1 

for the nth element of x (s^/* 

Linoaritv s ih terms of comoatable numbers, linearity of a 
digital signal processing algorithm would require that 
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(1) ct K(s^) x (s) = K(s^) a X(s) 

(2) a X^s) +x 2 (s) = a x^s) + a x^s) 

where x^(s) and X 2 (s) are sequences of input and K(s^) is a 
computable number. 

* Digital signal processing system is causal if the 
computable input values are as follows. 

x^n^s) = x^n^) for n <k 
then y-|_ (n,s) = y 2 (n^s^) for n < k 

even for x^Cn^s) & x 2 (n / s 1 ) for n > k. 

Here x^(n,s) and x 2 <n,s) ore input valuesond y-^Cn^s) and y 2 (n / s 1 ) 
are the corresponding output values. 

Stability : A digital signal processing algorithm is said to 
bo s table if any bound input values results in a bounded output 
values i.e. 

if x(n/s) < M x ( s i) 

for all values of n, then there is a computable number M y (s 2 > 
such that 

y(n,s) < 

for all values of n. Here M^s^ls a computable number 
which is greater than input values for all n. 
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6*3 Finite, word L.e^g^th a.nd Computable num bers 


Let us now examine the limitations resulting from 
finite word length from the point of view of computable numbers. 
Finite word length in the implementation of a digital hardware 
can produce three different type of errors. For the sake of 
eoncrctonoss let us take a specific case of a digital filter 
of first order which can be represented in real analysis as 
a difference equation. 


y(n) = x(n) - a.y(n-l) ... 

where y(n) and x(n) represents outputs and inputs of the digital 
filter. The difference equation when written in terms of 


computable numbers is as follows - 

y (n, s) = x(n,s 1 ) -a <s 2 ) . y(n-l),s 2 ) ( 6 * 2 ) 

The three different types of errors are as follows 

1. A real number x(n) when represented on a machine will have 
the so called input quantization error. 

2. The multiplier value <a‘, will be represented with the 
multiplier coefficient quantization error. 

3. in the rounding operation in multiplier, there exists 


rounding error. 

. . , . err or bound on the computable 

Total error wil3 be the error 

v . nnm nutable number concept also & cates that for 

number y(n,s), Cortputaox 
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value of ‘s’ once specified, there exists a method to compute 

Y (h/S> • For this we require error computation on y(n, s) and 

controlling it to a value less than K Control of each of 

s 

three error contributions can be achieved as follows : 

1 . To reduce input quantization error, decrease the step of 
quantization. 

2. To reduce multiplier coefficient quantization error, one can 

represent multipliers with greater precision, 

* 

3. To reduce rounding erroj? one can use longer word length 

For the same of simplicity we assume that each type of error 
con bo studied independently of the others. 

6.3.1 Input q uantiz atio n error s 

Lot us take A s — p as the step size for the quanti- 
zation o £ input and digital filter represented by difference 
equation (6*2) starts operation at n = 0 with initial condi- 
tions as follows. 

x (-1 ) « x(-2) = ... = 0 
y(-l) = y(-2) = > * » • “0 

So equation 

y(n,s) = x(n,s x ) - a(s 2 ) .Y (n-*l,s 2 ) 
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can be written for n » o as 

y(°/S 1 ) = x(0,s 2 ) .. a .O 

(since a(s 2 ) - a and no rounding error assumed) 


where 


For n = 1 


y(l/S 3 ) - r(l,s 2 ) - a.y(0 / s 1 ) 


where 


(1+a) 2" 


For n - k 


y(/S. i) x(k / s 2 >. a.yOc-lfSfc) 


whore 


frU 

s 2 S 


(1 + a + ••• + a k "" 1 2 ) 


2 , K-x. 

(1 -rcl'HSL + • • * & ) 


1 

for a <1 


( 6 . 3 ) 
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Thus the error contribution in the output at n=k 
is given by equation ^ . This indicates that after 

starting operation at n = 0, error is accumulating but once k 


is large,, a 11 << 1 as a is less than unity, in general* 


k+1 


1 

rr a ' 


Now let us suppose at n = k+1 input is removed* 

yCk-fl, s k+2 ) = - a *y(k, s^) 


where 


1 

■ ■i *■. ■ x 

^c+2 


a. 


°k 

1 a* 


Also 


Y y(k+2, s k+3 ) * -a.y (k+1 , s k+2 ) 


where 


s k+3 


2 p / l««a v 

= a • 2 ( Tzr* 


Thus 


y (k+m, s k+ i +rn 


) « -a*y (k+m-1, s^) 


m 


=* a 


^k-Hl 'Hn 


k 

, l~a ^ 


where 
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For a < 1 t 


is less than 


s' 


and 


is 


3 k+2 k+1 k-:-l nm 

stall lesser. This means that the quantization error contri- 


bution decays after the input sarrples are made zero. 


6 , 3.3 Mul tipl ier coefficient Quantiz^ .error jand Sensiluvit^: 

This error is contributed by approximate representation 
of a multiplier. Errors due to rounding in multiplier and input 
quantization is taken to be zero. Starting computation of 
y(n,s) using equation 6.2. * 

At n = 0 , initial conditions are as follows. 
x(~l) = x(- 2 ) = x(- 3 ) = ... - 0 
y (- 1 ) = y(-2) = y (- 3 ) = ... =0 

Then equation 6» 2 at n = 0 reduces to the form 
y (0, s 0 ) “ x{0) — a (S2) .0 

where — = 0 ; (x(0,s) - x(0) as assumed) 

s o 


At n -1 


where 


(1,82) « x(l) - a(s 2 ). y(0^ s Q ) 

1 1 ( 

—i— » max | a(s 2 X .y (°/ S D ) *-(a(s 2 )t . (y o / s 0 ^i s 


Y(o, s d ) 


x(0) 

Sn 
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At n> ~2 


y(2,s 4 ) * x(2) - a(s 2 ) .y(l/S 3 ) 


where 


lats,) .y(l,s,) ~(a(s 2 ) ) , (y(l,s,)+ ~ ) S 


3'~ s. 


1 / 1 v a ^ s 2' > 1 

?2‘ y< '"s’ + ” 3 ~ + "~®2* s 3 


x(l) - a(s 2 ) «x(0) 


x(0) 


x(0) 


+ a(s 2 ) .~— + T" 


x(l) x(0) 2 x(0}a(s 2 ) 


s 2 *" s 2 


(taking maximum value) 


5iil + x(0) ( (a(s 2 )-h ~ ) 2 -a 2 (s 2 )) 

^ s.-, 


For n = k 


Y (k# s^_^ 2 ) ~ x(k) — a(s 2 ) ,y(k— 


where 


3 k+2 m=0 


, \ t < \ , 1 \k-m k-m, v 

x(m) (a(s 2 )-i- ~ ) - a (s 2 ) 


« (6.4) 


This indicates that output error increases as the value 
of n increases. Another form of difference equation can be taken 
as follows . 

y(n,s) b 1 (s 3 ).x(n / s 1 )-a(s 2 )-y(n^l / s 2 ) 
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In this case now we can see output error contributed 

by error ~ in representation of b 1 (s.,) . Initial conditions 
s 3 1 J 

are given at as follows. 

x(~l) = x(-2) - . . . = 0 
y(~l/s) = y (-2^ s) = ... = 0 

Here input quantization erro rounding erro and error due to 

1 

— with a(s 0 ) is taken to be absent. 
s 2 2 


At n « 0 


y(o,s Q ) * b 1 (s 3 )x(o) - a.;jr{-i,s 2 ) 
= b 1 (s 3 ).x(0) 


(Here — =0 assumed) 


where 


o 


I 


At n 


y(l / s 4 ) = b 1 (s 3 ) .x(l)-a.y(o / s 0 ). 

= b 1 (s 3 ).x(l)-a.b 1 (s 3 ).x(0) 


where 


1 

s, 


x(l X 


ja 

s. 


X(J1 + 
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At n = 2 

y(2,s 5 ) 

where 

1 __ 
s 5 

Thus for n - k 

y(k, s k+3 ) = b i ( s 3 > *x(k)-a*y0o-l / s k+2 ) 


= ( s ^ X-xX 2 } i— a.y(l,s^) 


*S 2 l 


a 

+ =• 
&A 


x( 2 ) ^ 


+ 


^•*{21 
s.. 


where 


JL 

s k+3 


~ (x(k~l)+a.x(k-2) + 

s 3 



, v k~m 

£ x(m)a 
n-0 


-i- a^.xCo) ) 


(6.5) 


Here too, the error contribution increases with increase in n. 
This error can be reduced by reducing error bound 1/ s^ on 
b^(s 3 ) and l/s 2 on a(s 2 )» 

Sensitivity 5 A > concept related to this error contribution is 
the sensitivity of a digital filter structure. Sensitivity of 
a structure is defined as the ratio of error contribution in 
output due to change in mutliplier coefficient to change in 
multiplier coefficient. Structures taken are direct, cascade 
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and parallel form and each stage is restricted to first order 
only %/hose difference equation representation is given by 


equation 6*2 # 

1, Direct Form : Let a direct 
in Fig. 6.1(a) be 

y(n,s 1 ) - a c (s q1 ) x ( n ) 


form of a digital filter as shown 

aiCs^j) x(n-l) + ... + a^(s Nl ) 
.X(n-N) ... ( 6 .6 ) 


Let us take a (s -) multiplier value that have — 

O Ol 2 * * S ol 

error bound. All other multipliers are assumed to have no 
multiplier quantization error. Thus we can rev/rite equation 
( 6 . 6 ). 

y (n, s-j_ = a 0 (s ol )x(n)-,'-a 1 .x(n-l) ... + a^xCn-N) 


where 

Thus 


_1_ _ _ xCn)_ 

S 1 s ol 

1M 

sensitivity - - g- r ~" " 

s ol 


x(n) 


( 6.7) 


2. Cas cade Form s This form of structure is shown on figure 

6.1(b). All the stages are put in cascade. Let us represent 
th 

i stage as follows 

y^(n,s) — b^(s^)x(n) ~a^(s 2 )y(n~l*s 3 ) . 
where i = 1,2, ...N. 
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At first let b ± (s 1 ) have 1 /s. error bound and l/s_ = 


2 = 0 on a a(s 2 > 


Here 


y i (s) =x i+ i (s > 


or a Yitnjs) - = | ( xpnjsHa.. (s,)x, (n_l,s) 


x 2 ' i 

+ ... + a?(s 2 )x i (o <r s)) (6.8) 


from equation ( 6,5 ) 


This is now error in input (n, s) being denoted as E^, Thus 

output error 


Ay i+1 (n,s) = b i+1 ( E n + a 


x-KL n -1 


+ 


n-1 

a i+l *1 


This is again an input error to i+ 2 nd stage and so on. 
This may any error from first stage gets enlarged by successive 
stages, in other worda front stages are more sensitive t han the 
latter ones. 


Now for second case, let a. (s„) have + - error bound 

i / " s 2 

but l/s-^ = 0 on b^(s^) , Then using equation 5^4 we have 

n-1 1 

Ay i (n,s)-S n ^ S x n (s 1 a ± (s 2 )+ | ) nwIn -aj" rn (s 2 ) 

m«o z x 

* ♦ * ( 6 * 9 ) 


Hare y i (n,l) = x i4l (n,s) 
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or x i-KL^ rr,s ^ ^ las ± 3 n as e ^or/ which will result 

(n,s) 


AY 


i-rl 


(n,s) » b i+1 ( 6l > C 


E n + 


a i+l (s 2^ 


3 - 4 * 

n-1 


+ a' 


n~l 

ill 


(s 2 ) 


E l) 


Thus if a multiplier coefficient varies in first stage, output 
variation will be larger than a case when a multiplier coeffi- 
cient vary in a latter stage. 


3 . Parallel Form : This form is shown in Fig. 6 . 1 (c). All stages 
are put in paralle. Let us represent tms structure as follows 

N 

y (n, s) « E y.(n,s-> 
i=l 


and y i (n # s 3 ) = In (s 1 )x(n / s 1 > -a i (s 2 )y i (n-l,s 2 ) 


If we take the case when l/s 1 error on b ± (s^) exist alone. 

Ay(n,s) =^y^(n,s 3 ) = ~ ( x(n, s) -t-a^(s2) x(n-l ,s) +. . .+a!? (S2) x(o, s)) 


from equation ( 6« 5 ) 

1 


Simila rly if we have only a, (s 2 ) with error ~ , 

i ^ 2 

n “ 1 , . , , v . 1 v n-nr n-m 

AY (n, s) =Ay i (n,s 3 ) 


. . . . < X Il-UI, V 

2 x(m,s) (a i (s 2 )+ ~ ) -a ± (s 2 ) 

2 


irr =o 


from equation ( 5,4 ) 
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Thus here cascading effect is not there as was in cascade form. 
This may cascade form contributes more error at output due to 
a variation in multiplier coefficient of any stage except la st 
one, than the parallel form. Direct form has minimum sensiti- 
vity in these three forms. Parallel form is less sensitive than 
cascade form. 

6.3.3 Round inq _ _e_r_ror_ of Mu ltipl iers_ and Dea d band aff ect 

This error contributed in output is by multiplication 
of two numbers represented in b^ bits each rounded to b^ bits 
only and not 2b^ bits. Filter is represented by equation (6-2) 

y(n,s) = x(n,s 1 ) - a(s 2 ) .y(n-l, s 2 ) 

Here errors due to input quantization and multiplier coeffici- 
ent quantization are assumed to be absent. Say multiplication 
a(s 2 ) .y (n-1, s 2 ) has + p/2 error due to rounding. 

S tarting conditions at n— o are as follows 

x(-l) « x (— 2 ) = x(~3) = ... - 0 
y(-l) = y (-2) - y (-3) = ... = 0 

Then equation ( 5 , 2 ) at n = 0 reduces to the form 
y(o,s Q ) = x(o,s 1 ) - a(s 2 ).0 
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where |=0 ( i = o aS assumed) 

o ^ 

Here x(o # s 1 ) * x(0) and a(s 2 ) = & 

At n = 1 

y(l,s 3 ) = x(l) 3- a.y ( 0 /S Q ) 
where l/s 3 - 0 + 

At n = 2 

y(2,s 4 ) = x(2)-a.y(l/S 3 ) 

where l/s 4 = 0 + a/s 3 + p/ 2 
=s p/2 (1+a) 


At n = 3 

y(3,s 5 ) = x( 3 )-a.y( 2 / s 4 ) 

where l/e 5 = 0 + a/s 4 + p/2 

2 X 

— p/2 (1+a+a ) 


Thus at n = h 

y s ic+2 

where 1/ s x+2 


S = xOO - a.yCh-l/S^y) 
p/2 (1+a-r •••+& ) 


1 -a 

p/2 , 


X 
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For a large value of k and a <1 

i/s x+2 ^p/2. rhr 

1 

So error contributed at output will remain at p/ 2 , value 

after a large * value of n. 

Now we will see , if at n=k we remove input, then how 
this erro contribution is ‘going to change. 

At n - k+ 1 , innut is removed. 

y(k+l,s k+3 ) = - a.y(Jt /S]c+2 ) — 


where 



= p/2 +p/2 a. (1-i-a-i-. . .+ a ) 
— p/2 (1 -fat ••• + a ) 


At n ” k +2 

y(k+2, s k+4 ) = - a.yCk, 5^3) 


where 



— p/2 t a/ 

2 k 

= p/2 (1-l-a-i-a +-.*+a ) 


Thus at n s k-l-m 

y(k+m,s^. hnH _ 2 ) = " a.y (k+m-1, s leHfn+ i) 
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where — — — ■ 

S k-HrH-2 


p/2 + 


_a 

S k-Hn+1 


= p/2 

“ P/2 


(l+a+ , 
l-a n ' ! ™ 


1 —•a. 


r 


k+m~l 

a 


) 


So for n taken larger than 10 and a <1 


k-t-m+2 


f p/2 


1 

1-a 


( 6 . 12 ) 


Thus 1/ Cs. ^/ s ]t-f. m +? and this error is maintained though 
output will decay as input has been removed. By equation ( ) 

we see that yCk/S^..^ ■‘■ s b y a value less than 1 so 

we can write 

y(^/Sk +2 ) >y(k+1 ' s k-r3 } > y(k+2/ s k+4 ) etG * 


In Digital Signal Processing theory/ this is termed 
as the Dead band effect/ that is even if input is zero, output 
will oscillate between between p/2 . 1/1— a to -p/2.1/l-a. Let 
us now compute dead band effect bound using the conventional 
real analysis. 


The Dead-Band ef fect : We generally represent a digital filter 
by a difference equation of the form* in real analysis as follows 

n. n-1 

y (n) - S b.x(k) - E a-.y(l) 
k=0 iv 1=0 x 
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Let, us hove following conditions ; 

1. initially x(k) and y(l) are zero so y(n) = o. 

2. At n=0 input was non zero and one gets output also non zero. 

3. After a long time input is removed at n > n^. 

In such a situation the output y(n) for n > 1 ^ will 
also decay toward zero if system is stables and real analysis 
is used. 

In practice due to finite word length of the machine 
one will have rounding errors accompanying the output. 

For example for first order difference equation 

y(n) = x(n) - 0.8 yin-l) 
find the nature of y(n) 7 n > 0 

Given y(~l) = 8 , x(k) » 0 for k > 0 / Result in tabular form is 


n 

x(n) 

y(n-i) 

-ay (n-1) 

y(m) 

0 

0 ' 

8 

-6.4 

-6 

1 

0 

—6 

-!-4 . 8 

5 

2 

0 

5 

-4.0 

-4 

3 

0 

-4 

3.2 

3 

4 

0 

3 

-2.4 

-2 

5 

0 

-2 

1.6 

2 

S 

0 

2 

-1.6 

-2 



Here we carryout rounding to the n earest integer at each 
computation. The Dead band effect is visible after n=4 in 
the table given. 
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A general technique to study the dead band effect 
was given by Jackson [10,P296] Let us briefly discuss this. 

It is assumed that at each product a.y(n-l) in 
y(n) = x(n) - a.y(n-l) .. (6.1^ 

is rounded to the nearest integer. According to this rule 

Qnt. [ja.y (n~l)| ] = Integer [ ja.y(n-l)! -HD. 5 ] .. (6.1 4 ) 

When input is absent, 

y(n) = -a.y(n-l) 

The condition for y(n) and y(n~l) to be equal is that 

Qnt. [ i a.y (n-l)i ] = jy(n)| 

« |y(n-l)| 

or Iftt. [ |a.y(n~l)| + .5 ]= y(n-l) .. (6.15) 

from rule (6.14 ) 

or Jut. [ |y(n-l)| ~(i- la| ) iy(n-l)I + 0.5 ] =|y (n-1) I 

This means that 

0 < - (1- 1 a ! ) |y (n-l)l + 0.5 <1 
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or | y(n-l) j < .. (6.16) 

As y(n~l) can have only integer values for 
j a j < 0,5 ; y(n-l) = 0, Thus there is no dead band effect. 

Comparing equation ( 6.12 ) and equation ( 6.16) 
we find that computable analysis provide an explanation of the 
dead band effect. Now let us see whether the absence of dead 
band effect for a < .5 can also be explained using the concept 
of computable analysis. Answer is that for a < 0.5, the 
rounding error will not be there in such cases, thus so round- 
ing error will also decay as the value of y(n,s)« 



Section 7 


CpmiuAiPji 

The aim of this the *4.* has been to present the idea 
of computable number and the basic concept of computable 
analysis to those who are concerned about DSP. Originally 
it had been envisaged that whatever the results on computable 
analysis were already available/ they could be used to reformu- 
late the theory of digital signal processing in a straight 
forward manner. It however turned out that the task was not 
as straight forward and that a considerable amount of back- 
ground needed to develop, before it could be seen how compu- 
table analysis could be incorporated into Digital Signal 
Processing. The results of the efforts made in this direction 
are what is contained in this thesis* 

Among the major points that seem to emerge from- our 
study are the following 2 

1# A more satisfactory theory of Digital Signal Pro- 

cessing can be developed if we abandon real analysis. 

We use a kind of theory of numerical analysis that makes 
use of the theory of recursive functions and idealized 
computer model such as URM. 
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2 , in practice a lot of Digital Signal Processing is carried 
out using dedicated, hardware that has severe limitation 
of word length. The computable number analysis in this 
context provides a more appropriate mathematical framework 
within which signal processing algorithms can be studied. 

3. For the theory of computable analysis to become practica- 
lly usable, a great deed of work needs to be done to 
rewrite library functions and procedures and subroutine s 
for arithematic operations and relations in the language 
of computable analysis. 

Special direction for future work as the followings s 

1. Complete ore formulation of the theory of digital signal 
processing can be attempted, for which a start has been 
made in this thesis. 

2, Once we accept the idea of computable number as® a programme, 
an explicit algebra of programmes needs to be evolved. In 
this algebra given a programme we should talk of inverse 
programme and soim. The purpose of such an algebra will 

be to simplify the task of inplem anting computable analysis. 
Thus a computable number a(s) is given in form of a programme 
we should be able to cuickly construct the inverse programme 
corresponding to l/a(s). 
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GETTING S3UARE ROOT WITH 1/N ACCURACY 

BY 0 pSfSSfl T SlftiH*pSgn?Si A lg? 8 gi , SIIi IN ANALYSTS 

DIMENSION X(5,5),Y(5,5) ,Z(5) 

OPEN (IJNlTrS, r>EVICE='’ D5K * , FILEs'pmjT .FOR ' ) 
2P^NCUNIT=4,DEVICE='05K% FTLE= 'PIN. FOR') 
RE4DC4,*)NTRV,ERRNn ' 

CALL ADSU8(NTRM f ERBNn.ERBNDl ) 

ZTsQ, 0 

WRITE (5*2002) 

FORMAT(Sx,'Ya,l)',iOX,'YCI,2)',5X,'Ya,l)*YCI,2)',10X,'Si 

1 # 5 / 3 
TNBPsO 

DO 100 Isl,NTRW 

READ(4,*)CX(I,J),J=l f 4) 

^ i *** s, 0 

IF(i,0/S2«*ERBND 1)121,121,101 
S1=(Y(I,1)+¥CI,2) ) /ERBND1 

IF (|l«2 §8 00*0)1 02, 102,121 
7,1 sZT+ZCl) 

«RI|EC5^»)Y(I,l),Ya,2),ZCI),Sl,S2 
VRTTE (5 72005 1 ZT.ERBND.NTRM. inbr 

F0RHATC10X, 'TOTAL VALUE' . Fi 2. 3 , / , 10X , 'ERR, BOUND' ,F1 2 . 3 , / , 
1 ERMS ' » 15 , / , 1 OX , 'FEEDBACK LOOP , 16) 

STOP 

END 

SUBROUTINE TO CALCULATE £X)**1/N 

*0=0UTPUT AlsSTEP OF ROUNDING (2*S) AB=POWER( VALUE OF N) 
SUBROUTINE SQRTTCAO, AVALU,AB,A1 ) 

A2=0.0 

A3=A2/A1 

A4=A3**A8 

A5*A4-AVALU 

IF£A5)16,16,15 

A0=A3 

RETURN 

A2«A2+1.0 

GO TO 11 

END 

SUBROUTINE MULTIPLY El (BN)sXi 1 ( AN) *X22 ( AN) 

ONLY TO CALCULATE ERROR BOUND "BN* 

SUBROUTINE MLTPYCX1 1 ,X22 , AN , BN , El ) 

FF=1.0/AN 
ElsXi 14X22 
FRTTsFF* CX11+X22) 

BN=1 ,0/ERTT 
RETURN 

JTMf) 

SUBROUTINE DIVIDE FICC»)«X33CCJf)/X44CCN) 

ONLY TO CALCULATE ERROR BOUND *DN* 

SUBROUTINE 0VDE(X33,X44,CN,DN,FI) 

GF=1,0/CN 

FI*X33/X44 

EROR*GF4C1,0/X33+1.0/X44) 

DNsi.O/EROR 

RETURN 

SUBROUTINE ADDITION/SUBTRACTION HI (EN)s X 5 5 (EK) +/-X66 (EK)+/ 

ONLY TO CALCULATE ERROR BOUND "EK" 

SUBROUTINE ADSUBf N,EN,EK) 

EKSEN/FLOATCN) 

RETURN 

CALCULATION OF SINCX(S)) WHICH HAS +/-(1/S3 ERROR BOUND 
BY PRITHVI SINGH PUROHIT 
SUBROUTINE SINE(X,5l ,S2,SNE) 

ENsl.0/32 

SNEsX 

Nsl 

SNXsSIN(X) 

T2*X 

T3»X+EN 

SIXE*SIN(T3) 

SNERsX+EN 

NsN+1 

SNE*SNE+T2 

SNERsSNER+T3 

EN1sT2*R1/(1. 0"R1? 


',10X 


10X,'T 


uu 


205 

201 

2001 

200 «? 

200 


ETT=ABSCEN1}+ABS(EN2) 

IF CABS (EMI ) -ABSCEN2) >201,205, 205 

STTaJ ,0/ETT 

IF(S2*STT)200,200#201 

IFCN-200001 202 ,202,200 

WRITE (5, 2001) SNE,SNER,N,ETT,ENt,EN2 

WRITE(5,2000)S?JX,SIXE,X,S1,S2,STT 

FORM AT f2F12. 6, 112,^12.7) 

F0FMATC2F12,6,4F12.3,//) 

RETURN 

END 
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SSXSiYI SINGH PU R OHIf EE/8310421 
PROGRAM TO OBTAIN INVERSE OF A MATRIX 

t-S ▼ K"' »T <S» 'T n ®LT n/4 A 4 a * %#**>-%*.** - . V V- * A *> 


30 

101© 

42 

41 

44 

31 

32 


34 


35 


38 
33 

C 4« 
40 

39 


36 

3003 

37 

3000 

88 



ER C 10 , 203 

, JE.VICKS'nSR'.FTT.Ii’s'MTM MD'I 

READ ( 4 , *) N ,S 
51*1,0 
NCNTRsO 
DO 30 1 = 1, N 

READC4,*)CB(I,J),J*1,N) 

WRTTEC $ #1010) 

FORM AT C IDX, 'MATRIX "A"') 

DO 42 1 = 1. N 

WRITF(5,*)CB(X,J),J=1 ,N) 

Mssfl+N 

W?=N +1 

NCNTR=NCNTR+1 
DO 31 LI*1,N 
DO 44 LC=i,N 
A(LT f LC)=B(LI,LC) 

DO 31 LJ=M2,* 

FRCLI,LJ)=0,0 
AC t*I ,1.4 3 *0,0 
no 32 K* 1 ,N 
I2*K+N 
ACK, 123*1.0 
ERR=0.0 
DO 33 LJ=l , M 
•J2=LJ + 1 
PsACltJ.IiJ) 

DO 34 1 = 1 , M 

ERCLJ,I)=CABSCACLJ,I))+1.0/Si)/CABSCP)-1.0/Sl)-ABSCACLJ,I)/P) 
A( LJ , I 3 =A (LJ, I) /P 
DO 33 liKsl.N 
DO 33 LI=J2,M 
IFCLK-LJ 3 35,33,35 

ERCLK,LI)=Cl*0+ABS(A(LK,LJ) )+ABS(A(LU,Ll) ) )/51 
ACLK,LT)*ACLK,LI)-ACLJ,LI)*ACLK,LJ) 

IF (ERR-ERCLK, LI) 338,33,33 

ERR=ER(LK,LI J 

CONTINUE 

IF C ERR" 1 • 0/S) 39, 39,40 
S1=S*S1*ERR 
5 1 *S 1 + 1,0 

IFCSl-10000.0341,41,39 
DO 36 1*1, N 
DO 36 J=M2 , M 
L3*J-N 

ERCI ,L33=ER(I» J) 

AI(I,L3)=A(I. J) 

FORMaIc^j/?10X, 'INVERSE "A" MATRIX') 

FORMAT C§X'VALUfe S 6 FERR 0 R IN ELEMENTS OF MATRIX "A* ' ,F 10 . 2 , 
l°5X,'MAxi, ERROR IN ELEMENTS OF MATRIX ■ INVE RSE A , 
1 5X, 'FEEDBACK LOOPS', 16,/,/, 10X, 'ERROR MATRIX OF INVERSE A 
DO 88 1 * 1 , N 

WRITEC5, * J (ERC I, J) , J*1 r N) 

STOP 

END 
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